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Abstract 



Let X be a smooth manifold of dimension 1 + n endowed with a 
qh Lorentzian metric g. The energy tensor of a 2-form F is locally defined as 

^ ; T ab := - (FjFu - i F^F ij9ab ) . 

bi) In this paper we characterize this tensor as the only 2-covariant natural 

tensor associated to a Lorentzian metric and a 2-form that is independent 
of the unit of scale and satisfies certain condition on its divergence. This 
characterization is motivated on physical grounds, and can be used to 
justify the Einstein-Maxwell field equations. 

More generally, we characterize in a similar manner the energy tensor 
associated to a differential form of arbitrary order k. 

Finally, we develop a generalized theory of electromagnetism where 
charged particles are not punctual, but of an arbitrary fixed dimension p. 
In this theory, the electromagnetic field F is a differential form of order 
2 + p and its electromagnetic energy tensor is precisely the energy tensor 
associated to F. 
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Introduction 



Let (X, g) be a relativistic spacetime of dimension 1 + n; that is, X is a smooth 
manifold of dimension 1+n and g is a Lorentzian metric of signature (+,—,.?.,—). 

An electromagnetic field on X is represented by a differential 2-form F, and 
its electromagnetic energy tensor T is a 2-covariant tensor defined in a local chart 
by the formula: 



The main purpose of this paper is to prove the following characterization of 
this energy tensor: 

Theorem: The energy tensor is the only 2-covariant tensor T = T(g, F) natu- 
rally associated to a Lorentzian metric g and a 2-form F, satisfying the following 
properties: 

1. T is independent of the unit of scale; that is, T(X 2 g, XF) = T(g,F) for 
any A > 0. 

2. At any point, F x = =>- T x = 0. 

3. If dF = then divT = -i dF F. 

More generally, the energy tensor can be defined for differential fc-forms ( [13] ) , 
and we characterize these tensors in a similar manner (Theorem 13.51) . 

Let us briefly explain the role of the above theorem in order to physically 
motivate the definition of the energy tensor of an electromagnetic field. In Gen- 
eral Relativity, the matter content of spacetime X is represented by a symmetric 
2-covariant tensor T m (the matter stress-energy tensor) and, in absence of elec- 
tric charges, the mass-energy and impulse conservation laws are encoded in the 
equation: 



Nevertheless, when dealing with charged matter, the Lorentz force law im- 
poses: 

divT m = ijF = i dF F 

where J is the charge- current vector field and dF = J* because of the second 
Maxwell equation (in the interior product, dF stands for the vector metrically 
equivalent; see our notations in Section II. ip . Therefore, in order to have the 
aforementioned conservation laws, it is necessary to assume that, apart from the 
stress-energy tensor T m of the matter distribution, there also exists some kind 
of energy associated to the electromagnetic field itself, represented by some 2- 
covariant tensor T e \ m , such that: 




divT T 



m 



. 



div (T m + T clm ) = . 
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Of course, this equality implies divT e i m = —igpF. Since dF = (first 
Maxwell equation), the tensor T c \ m has to satisfy condition 3 of the Theorem. 

As concerns to the first hypothesis, observe that if the metric g is changed by 
a proportional one, ~g = X 2 g ( A G K + ), then the proper time of the trajectory of 
any particle is multiplied by the factor A. In other words, replacing the metric 
g by ~g = X 2 g amounts to a change in the time unit. This change modifies the 
other units of length, mass and charge, since we assume they are chosen in such 
a way that the universal constants (light velocity, gravitational constant and 
Coulomb's constant) are all equal to 1. It is easy to check that a change of the 
time unit ~g = \ 2 g (with the corresponding change in the other units) implies 
a modification of the type F = XF in the mathematical representation of the 
electromagnetic field, while the matter tensor remains invariable: T m = T m . This 
last equality and the equations div (T m + T e i m ) = = div (T m + T clm ) , imply that 
the electromagnetic energy tensor also stands invariable: T e \ m = T e \ m . That is to 
say, the tensor T c \ m has to satisfy condition 1 of the Theorem. 

Finally, condition 2 of the Theorem states that the electromagnetic energy is 
null wherever the field is null. 

Summing up, these are three properties that have to be satisfied by any physi- 
cally reasonable definition of electromagnetic energy tensor, and our result proves 
that the choice is then uniquely determined. 

This problem of characterizing the electromagnetic energy tensor is classical 
and has already been studied in the literature ([I], [6], [7], [9], [ID]). The closest 
result to our statement is Kerrighan's ([I]), where the tensor T(g,F) is assumed 
to be symmetric and its coefficients are assumed to be functions of the coefficients 
of g and F (so tensors using higher derivatives of g and F are not considered). 

Both restrictions are removed in our theorem where, instead, we require inde- 
pendence of the unit of scale. This is a physically meaningful condition which, in 
spite of its innocent appearance, turns out to be very restrictive. Some examples 
may illustrate this point: 

- The Levi-Civita connection is the only linear connection naturally associated 
to a metric that is independent of the unit of scale (Epstein, [3]). 

- The Einstein tensor is the only, up to constant factors, 2-covariant natural 
tensor associated to a metric that is divergence-free and independent of the unit 
of scale ( Navarro- Sancho, [12]). 

- The Pontryagin forms are the only differential forms naturally associated to 
a metric that are independent of the unit of scale (Gilkey, [1]). 

Let us summarize the content of the article. 

We begin with a preliminary section where we recall the definition of energy 
tensor associated to a differential fc-form u and the main properties that we use 
to characterize it. 

In the following section, the problem of computing natural tensors associated 
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to a Lorentzian metric and a /c-form uj, subject to a certain homogeneity condition, 
is reduced to a problem of representations of the orthogonal group. The main 
result (Theorem 12. 141) is similar to other known results (|15j. |16j . [T7]). but with 
a slightly different homogeneity condition. 

Next, we determine all the 2-covariant natural tensors T(g,u)) that are inde- 
pendent of the unit of scale. As a consequence, it follows the announced charac- 
terization of the energy tensors (Theorem 13.51) . 

Finally, the existence of an energy tensor associated to a differential form 
of arbitrary order suggests the question of a possible physical interpretation for 
it. In the last section, we consider a generalized theory of electromagnetism 
for charged p-branes, introduced by Henneaux and Teitelboim ([5]), where the 
electromagnetic field F is a differential form of order 2 + p. We extend this 
theory up to the point of including fluids of charged p-branes; the corresponding 
Maxwell-Einstein equations require an electromagnetic energy tensor, which turns 
out to be the energy tensor associated to the form F. 



1 Preliminaries 

Throughout the paper, let (X, g) be a Lorentzian manifold of dimension 1 + n, 
whose metric has signature (+, — , — ). We assume X is oriented, with volume 
form dX, and time oriented. 

1.1 Notations and conventions 

Given a g-vector D\ A . . . A D q and a differential fc-form u, with q < k, we write: 
zdiA...ad,w := i-Dq ■ ..i Dl u = aj(D h D q , _) . 

Analogously, if u; is a /c-form and a) is a g-form, with q < k, we write: 



where oo* is the q- vector metrically equivalent to uj. 

With these notations, the metric induced on the bundle of /c-forms is: 

(u,u) := i^u = ■£U 3l '" 3h u jl ..j h . 
The Hodge star is the linear isomorphism *: A k T*X — > A 1+n_fc T*X defined 

as: 

*u := i^dX 
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and, with these conventions, it holds: * * cu = (— l)( fc+1 ) n a;. 

The co differential d: Q k (X) — > is the following differential operator: 

d := * d * , or, equivalents, * d := (-l) fc d * . 

In a local chart: 

(duj) il ...i k _ 1 = -V a u ail ... ik _ 1 . 

Remark 1.1. Later, we will need the following formulae for the components of the 
1-forms idujCJ and i^du: 



(k-l)\ 

#^- lfcV ^- + (^l)! 



{idwV)b = 77 —^aU al2 - ik Ubi 2 ... lk , 



{i u du)) b =^ 1 -^u il - ik V h u h ... ik + — — — cA^V^...^ . 



1.2 Energy tensor of a differential A;- form 

Let w be a differential /c-form on X. 

Definition 1.2. Let U be an observer at a point x (that is, U is a unitary timelike 
vector oriented to the future). Let us consider an orthonormal frame (-Do — 
U,Di, . . . , D n ) of T X X and the corresponding dual base (9 = U*, 9\, . . . , 6 n ). 

In terms of this basis, the /c-form uj decomposes as a multiple of 8 Q , called the 
electric part Ejj, and other terms without 8q, called the magnetic part Bjj\ 

co = Ejj + Bjj = (terms with 6 ) + (terms without 6$) • 

In other words: 

Ejj := U* Aiuu , Brj := iu(U* Acu) , 

so these /c-forms Ejj, Bjj depend on the observer U but not on the chosen basis. 
Moreover, as Ejj and B\j are orthogonal: 

(to, to) = (EutErf + iB^Bu) . 

These two addends have definite signs, that we modify to make them positive: 

^■^{-if-^Eu^u) = (-l^iiuuJuu) 

\Bu\ 2 := (-l) k (Bu,Bu) = (-l) k (U* A u, U* A a) . 

Hence, 

(w.w) = (-I)'- 1 (\Ev\ 2 - \Bv\ 2 ) 
and the right hand side of this equation does not depend on the observer. 
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Definition 1.3. The energy of a differential /c-form ui with respect to an observer 
U is the smooth function: 

e(C0:=i(|^| a + |B p | a ) . 
Unfolding the definitions, we see the energy is quadratic on U : 
e(U) = ~ (\E V \ 2 + \B V \ 2 ) = 1 -{\E V \ 2 + (-l) k (u,u) + \E V \ 2 ) 

= (W, it/w) - ^ (w, u) (U, U) 

so we are led to consider the corresponding symmetric tensor: 

Definition 1.4. The energy tensor of a differential /c-form u is the 2-covariant 
symmetric tensor T defined as: 

(-l)^ T {D 1 ,D 2 ) := (i Dl u,i D2 u)-^(cu,u)g{D h D 2 ) . 

This definition is made so that T(U, U) = e(U) for every observer U. In a 
local chart, 

(—l) k ~ l ( . . 1 . . 

ab ~ (jt _ 1)! I a W6l 2...* fe _ U H-ik9ab 

Remark 1.5. These energy tensors are a particular case of the superenergy tensors 
introduced by Senovilla ([13]): the superenergy tensor associated to a differential 
fc-form is precisely the energy tensor of Definition 11.41 

Next, we quote the main property of the energy tensors, although we will not 
use it in this paper: 

Theorem 1.6 ([H], Th. 4.1). The energy tensor T of a k-form satisfies the 
dominant energy condition. 

In other words, for any pair Ui, U 2 of observers (unitary timelike vector fields 
oriented to the future), it holds: 

T(U 1 ,U 2 )>0. 

For any observer U, the Hodge star maps the electric and magnetic parts of 
oj into the magnetic and electric parts (up to signs) of *u: 

*Eu(u) = ±Bu(*u) , *Bu(u) = ±Eu(*tu) . 

Therefore, u and *u have the same energy respect to any observer U and, conse- 
quently, both forms have the same energy tensor. 

The following three properties are easily obtained from the definitions, and 
they will suffice to characterize these energy tensors: 
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1. Let us write T(g,u) to indicate that the energy tensor depends on the 
metric g and on the fc-form u. For any A > 0, 

T(A 2 ^,A fc - 1 o;) = T(g,u) . 

2. At any point x G X and for any observer U, it holds e(U)(x) = if and 
only if = 0. Hence, 

T x = <(=> ^ = . 

Proposition 1.7. The energy tensor T of a k-form u satisfies: 

div T = i^du — ig^u . 
In particular, if u is closed and co-closed (du = = du), then divT = 0. 
Proof. In a chart, 

(divT), = V a T\ = V a ^■■■^.i fc -io; <1 -"V...* fc ^) 

= S ((Va^ l2 -^)^ 2 ... lfc + a;^- J HV a ^,.,J - i^^V^,., fc ) 

rm i ■ • f-i) fc ■ • 

^ -(iaww)t + (?u,du;) b . 

□ 

Corollary 1.8. For an?/ k-form u, the 1-covariant tensor (i_u,i_uj) satisfies: 

div = (— (^dw - ia u oj) + -d(a;,a;) . 

Proof. By definition of T, we have: 



Hence: 



div = (-l) fc ~ 1 divT+ ^ div ((w, w)a) 

= (-l) fc_1 (? w dw - ifou) + - d(w, u) . 



□ 



Remark 1.9. Analogously to the case of a 2-form, the energy tensor T of a A;- form 
a; also appears as the Euler-Lagrange tensor of a variational principle. Namely, 
if we fix a /c-form u; and consider the variational problem of order defined by 
the lagrangian density (u, u) dX on the bundle of Lorentzian metrics, then its 
Euler-Lagrange equations are precisely T = 0. 
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2 Natural tensors associated to a metric and a 
form 



In this section we study natural tensors associated to a Lorentzian metric g and a 
/c-form uj. The notion of "natural construction" is formalized within the language 
of natural bundles and natural operators. Our presentation slightly differs from 
the standard approach ([8]), so we give a brief exposition of it. 

2.1 Natural operators 

In the following, all bundles E — > X are assumed to be sub-bundles of some 
bundle of tensors on X (observe E — y X need not be a vector bundle; v.gr., the 
bundle of Lorentzian metrics). This hypothesis is not essential, but simplifies the 
exposition. 

Definition 2.1. Let E, F — > X be two bundles over X and let £, J 7 be their 
sheaves of smooth sections, respectively. A morphism of sheaves T : £ — > J 7 is 
called a regular operator if, for any smooth family {e s } se s of local sections of 
E depending on certain parameters, the family of sections {T(e s )} s£ s of F also 
depends smoothly on those parameters. 

According to a fundamental result due to Slovak ([H]), the regularity condi- 
tion of a operator T : £ — > J 7 implies the existence of a smooth map T : J°°E — > F 
such that T(e) = T o j°°e for all local section e of E, so that T(e) depends on 
the oo-jet of e. 

Definition 2.2. A bundle E — > X (sub-bundle of a bundle of tensors) is said to 
be natural if it is stable with respect to the action of local diffeomorphisms of 
X. 

That is, for any diffeomorphism r: U — > V between open sets of X, it holds: 

e e £(V) => r*e G £{U) . 

Definition 2.3. Let E, F — > X be natural bundles. A regular operator T: £ — > 
T is said to be natural if it is equivariant with respect to the action of local 
diffeomorphisms of X. 

That is, for any diffeomorphism r: U — > V between open sets of X, it holds: 

T(r*e) =r*(T(e)) . 
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2.2 Normal tensors 



When dealing with jets at a point, it is useful to write the "normal expressions" 
of the geometric objects under consideration. 

Let x G X be a point and let g be a germ of Lorentzian metric at x. 

Definition 2.4. A chart (zq, . . . , z n ) on a neighbourhood of x is said to be a 
normal system for g at the point x if the geodesies passing through x at t = 
are precisely the "straight lines" {z (t) = Aq£, • • • > 2Vt(t) = X n t}, where Aj G M. 



Remark 2.5. Via the exponential map exp x : T X X — > X, normal systems at x cor- 
respond bijectively to linear coordinates on T X X. Therefore, two normal systems 
at x differ on a linear transformation. 

Normal systems are characterized by the following well-known lemma 

Lemma 2.6 (Gauss Lemma). Let (zq, . . . , z n ) be germs of a chart centred at 
x G X . This chart is a normal system for a germ of a Lorentzian metric g if and 
only if the metric coefficients g^ satisfy the equations: 



9ij z j 9ij( x ) z j ■ 



j 



Definition 2.7. Given a a normal system (z , . . . , z n ) for g at x, let us write 

d r gij ( , 

9ij, ai ...a r ■— o p. \ x ) ■ 

u ^a\ ' ' ' u ^a r 

For any integer r > 0, the r-th normal tensor of g at x is defined to be 

9l '■= ^ 9ij, ai ...a r d x Zi ® d x zj <g> d x z ai <g> . . . ® d x z ar . 

Analogously, if w is a differential fc-form and we write: 

^i-\...iu,a\...a r ■ ~o ^ l-*v 

<9z ai . . . <9z as 

then, for any integer s > 0, we define the s-th normal tensor of a; at a; to be 

^2 u; h...i k ,a 1 ...a s d x z h <g> . . . <g> d x z ik <g> d x z ai . . . <g> d x 2: a3 . 



il,...ifc,ai,...a B 



Using the Gauss Lemma, it is easy to check that the normal tensors g r x have 
the symmetries stated in the following definition. 
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Definition 2.8. For each integer r > 1, the vector space of metric normal 
tensors of order r at x is the vector subspace N r C ® r+2 T*X of tensors P with 
the following symmetries: 

1. They are symmetric in the first two and the last r indices: 

P - P P - P w n • o 

1 ija\...a r 1 ]ia\...a r i 1 ija\...a r 1 U a tr(i) •••<W r .) v u *~ '-'r ■ 

2. The cyclic sum over the last r + 1 indices is zero: 

P + P + + P =0 

For r = 0, we define N := M x to be the space of Lorentzian metrics at x G X. 

Remark 2.9. Due to symmetries, N\ = 0, and therefore g\. — for any metric. 

Definition 2.10. For any integer s > 0, the vector space of /c-form normal 
tensors of order s at x is defined as: 

A s := A k T*X ® ^ S T X *X . 

Of course, for any fc-form a; we have 6 A s . 

Remark 2.11. Although we will not use this fact, let us remark that the se- 
quences of tensors (g x , g%, . . . , g r x ) and (g x , R x , (VR) X , {V r ~ 2 R) x ) mutually 
determine each other, as so happens with the sequences (cu x , uj x , . . . , uj x ) and (cu x , 
(Voj) x , . . . , (V s oj) x ), once the metric g is fixed. 

2.3 Computation of natural tensors associated to a metric 
and a form 

Let M C S 2 T*X be the bundle of Lorentzian metrics on X, let A k X be the 
bundle of differential A;-forms on X and let T«X = ® P T*X ® (g) 9 TX be the 
bundle of (p, g)-tensors on X. 

Their sheaves of smooth sections will be written, respectively, 

Metrics , Formst , Tensors^ . 

Definition 2.12. A natural (p, g)-tensor, associated to a Lorentzian metric and 
a differential fc-form, is a natural operator T : Metrics x Formst — > Tensors'^. 
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Definition 2.13. A natural (p, g)-tensor T is said to be homogeneous of 

weight weR if, for any metric g, any /c-form u and any real number A > 0, it 
holds: 

T(\ 2 g,\ k - 1 u) = \ w T(g,tu) . 

If T is homogeneous of weight 0, we say it is independent of the unit of 
scale. 

Let x G X be a point and let g x be a Lorentzian metric at x. We will write 
9x := 0(1, n) for the orthogonal group of (T x X,g x ). The symmetric powers 
S d N r and S C A S are linear representations of 9x . 

If V and W are linear representations of 9x , we denote Homo (V, W) the 
vector space of O g:E -equivariant linear maps V —> W. 

The following theorem allows to compute homogeneous natural tensors: 

Theorem 2.14. Let us fix a point x E X and a Lorentzian metric g x at x. There 
exists an M.-linear isomorphism: 

{Natural (p,q)-tensors homogenous of weight w} 

Hom 0te (S d * N 2 <g> • • • ® S dr N r <g> S C °A ® . . . ® S C *A S , (T^A) X ) 

where the summation is over all sequences of non-negative integers {d 2 , . . . , d r }, 
r > 2, and {cq, . . . , c s }, satisfying the equation: 

2d 2 + . . . + r d r + co + 2d + . . . + (s + l)c s = p-q-w . (2.3.1) 

// t/ws equation has no solutions, the above vector space is reduced to zero. 

This theorem is closely related to results of Stredder ([H], Theorem 2.5) and 
Slovak ([15], Theorem 3.3) and the proof is similar. 

Remark 2.15. If <p: S d2 N 2 ® ■ ■ • ® S^A/. ® S C °A <8> . . . (g> S Cfl A s (TjjX) x is 
an O^-equivariant linear map, then the corresponding natural tensor T(g,u) is 
obtained by the formula: 

T{g, u) x = ((^<g> .«fe g> ^) <g> . . . ® (£® (g) g r x ) <g> . . . ® ( w * <g> <g> u s x )) 

where (g x , g x , . . .) is the sequence of metric normal tensors of g at x and (o;°, ui x) . . .) 
is the sequence of fc-form normal tensors of u at x. In this equality, g is assumed 
to have the prefixed value at x. 
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Remark 2.16. The 9;c -equivariant linear maps that appear in the theorem can 
be explicitly computed using the isomorphism: 

Rom 0gx (S d *N 2 ® ■ • • ® S C *A S , (T*X) X ) 

Hom 0te (S d2 N 2 ® ■ ■ ■ ® S C °A S ® (TJX) X , R) . 

and applying the Main Theorem of the invariant theory for the orthogonal group 
Og x (see [2], Th. 4.1, for a proof in the Lorentzian case). This theorem states that 
any 5;c -equivariant linear map S d2 N 2 <E> • • ■ <8> (T£X) X — > R is a linear combination 
of iterated contractions with respect to the metric g x . 

As an example, for a non zero linear map to exist, the total order (covariant 
plus contravariant order) of the space of tensors S d2 N 2 ® • • • <8> (TgX) x has to be 
even. 

3 Characterization of the energy tensors 

In this section, we characterize the energy tensor of a fc-form by three conditions 
(Theorem 13.51) . To do so, we analyse separately the consequences of each of these 
conditions. 

Proposition 3.1. Let T : Metrics x Formsk — > Tensors® be a natural tensor, 
with k 7^ 1,3. 

If it is independent of the unit of scale, then T(g, u) is an H-linear combination 
of the following four tensors: 

Ricci(g) , r(g)g , , (w,w) g g 

where r(g) stands for the scalar curvature of the metric g. 

Proof. By Theorem 12.141 such an homogeneous natural tensor of weight w = 
is determined by a 5i -equivariant linear map: 

S d2 N 2 <g> . . . <g> S dr N r <g> S C «A <g> . . . <g> S C °A S ► (g) 2 T*X 

where the coefficients dj, Cj G N satisfy the equation: 

2d 2 + . . . + rd r + c + 2c l + . . . + (s + l)c s = 2 . 

If some Q is non zero, then there are only two possibilities: 

• Co = 2, ci = . . . = c s = dj — 0. In this case, we are reduced to compute 
9:E -equivariant linear maps: 
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S 2 (A k x X) ® T X X <g> T X X — ► R . 

As explained in Remark 12.161 those linear maps are linear combinations of 
iterated contractions. Due to symmetries, any such an iterated contraction is a 
linear combination of these two: 

T i y T ai akCll ak bb i T i y T ba2 akCa2 ak i) C . 

These contractions, in turn, correspond, respectively, with the tensors: 

(u,u) g g , (i-u,i-u) g . 

• ci = 1, Co = C2 . . . — c s = dj — 0: Any such a tensor would produce a 
9:r -equivariant linear map: 

A*X <g> T*X <g> T X X <g> T X X — )• R . 

but there are no such maps for k even or k > 5 (the contraction of two skew- 
symmetric indices is zero). 

Finally, if the q are all zero, then the tensor T(g,u) does not depend on u 
and therefore it is a linear combination of Ricci(g) and r(g)g (see details on 
[T2j . Theorem 5.1). 

□ 

Remark 3.2. In the previous Proposition, if k — 3 there also exists the 9a ,- 
invariant linear map: 

A 3 T*X <g> T*X <g> T X X <g> T^X — >• R , T ^ T abcabc 

that corresponds to the natural tensor Coi(Vw), where C i denotes the contrac- 
tion of the first two indices. 
If k = 1, 

t;x g) t;x ® t x a <g> t;x -> k 

there exist three different iterated contractions: 

that correspond, respectively, with the natural tensors: 

(divgoj)g , V 9 w , (VgUjf 
where (V^)*^, D 2 ) := (V g u)(D 2 , D x ). 
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Proposition 3.3. Let T: Metrics x Formsk — > Tensors® be a natural tensor. 
If it satisfies: 

1) T is independent of the unit of scale: T(\ 2 g, A fe_1 u;) = T(g, u) for all A > 0, 

2) At any point, u x = =>- T(g,u) x = 0, 

then there exist universal constants ^i^GM such that: 

T(g,u) = n 1 (i_u,i_uj) g + fi 2 (u,uj)gg ■ 

Proof. Condition (2) rules out the tensors Ricci(g) and r(g)g in the previous 
proposition, as well as the other exceptional tensors in the cases k = 1, 3. 

□ 

Theorem 3.4. If a natural tensor T : Metrics*. Forms k — > Tensors® satisfies: 

1. It is independent of the unit of scale: T(X 2 g, A fc_1 u;) = T(g, u) for all A > 0, 

2. At any point, uj x = T(g, u) x = 0, 

3. dxv g T(g, aj) = whenever u is closed and co-closed, 
then T(g,u) is a constant multiple of the energy tensor: 

E(g,u) := (-l)^ 1 ((i-u),i-u)g - ^(u,uj) g g 

Proof. By the previous proposition, there exist universal constants fi\,fi2 £ K 
such that: 

T(g,u) = ni{i-u,i-u) + H2{u,w) g . 
Then, writing T = T(g,u), 

divT = div ({i-co, i-w)) + /i 2 div ({cu, cu) g) 

/ii ^(-l)* 1 " 1 {i^du - i dul uj) + ^d(w, u}j + u, 2 d(uj, lu) 

By hypothesis (3), if u is closed and co-closed, then divT = 0. Comparing 
with the previous equation, we have 

= div T = Hi ^0 + ^ d(cu, cu^J + /i 2 d(cu, u) , 

hence \xi = — /xi/2 and we conclude: 



T = /xi(i_w,i_w) - — (cu,cu)g = t l i(-l) h 1 E . 



□ 
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We may reformulate the above theorem so as to eliminate the constant factor: 
Theorem 3.5. If a natural tensor T : Metrics*. Forms k — > Tensors® satisfies: 

1. It is independent of the unit of scale: T(X 2 g, X k ~~ 1 u) = T(g,u) for all X >0, 

2. At any point, u x = =>■ T(g,cu) x = 0, 

3. dw g T(g,u) = —iq^oj whenever uj is closed, 

then T(g,u) coincides with the energy tensor E(g,u). 

Proof. It is clear that T satisfies the hypotheses of the previous theorem, so T 
is a constant multiple of the energy tensor E. As both tensors have the same 
divergence whenever to is closed (see Proposition 11.71) . that constant has to be 
one. 

□ 

Remark 3.6. Let ClosedFormsk be the sheaf of closed fc-forms on X. We state, 
without proof, another variation of the previous result: 

Theorem. Let T : Metrics x ClosedFormsk — > Tensors® be a natural opera- 
tor. If it satisfies: 

1. It is independent of the unit of scale: T{X 2 g,X k ~ 1 u) = T(g,u) for all X > 0. 

2. At any point, u x = T(g,u) x = 0. 

3. div g T(g,uj) = -i duJ uj. 

then T(g,oo) coincides with the energy tensor E(g,u). 

4 Electromagnetism of £>-branes 

There exists a theory of electromagnetism for charged p-branes ([5]), where the 
electromagnetic field is represented by a differential (p + 2)-form F . In the rest 
of the paper, we shall extend this theory up to the point of including a force law 
for fluids of charged p-branes and an electromagnetic energy tensor, necessary to 
state the Einstein equation. This tensor is precisely the energy tensor of the form 
F introduced in Definition 11.41 

In this section we analyse the interaction of a charged p-brane with an ar- 
bitrary electromagnetic field. Our analysis is developed at the classical (non 
quantum) level and, in contrast to [5], it is based on the elementary concepts of 
impulse and acceleration of a p-brane. 

From now on, let us fix an integer p, such that < p < n, and let us write 
q := n — p. 
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Definition 4.1. The trajectory of a p-brane is, by definition, an oriented 
smooth submanifold S C X of dimension p + 1, whose metric g\$ has signature 



Associated to any p-brane, we also assume two constants, called tension t > 
(or mass, in the case p = of punctual particles), and electric charge q G M. 

4.1 Impulse form of a ]>brane 

In absence of external forces, the trajectory of a punctual particle is a geodesic of 
spacetime. To extend this fundamental principle to the movement of a p-brane, 
let us recall two different characterizations of geodesies: 

1. The trajectory of a particle is a geodesic if the impulse vector mil is parallel 
along the trajectory (where m is the mass of the particle and U is the unitary 
tangent vector to the trajectory). 

2. The trajectory of a particle is a geodesic if it minimizes the action m j dr, 
where r is the proper time of the trajectory. 

To determine the movement of a p-brane in absence of external forces, it is 
common in the literature to follow the second approach, using variational prin- 
ciples. To be precise, the generalized action is the Nambu-Goto action, t L dS, 
where dS is the (p + l)-volume of the trajectory S of the brane. 

Instead of that, in this paper we generalize the concept of impulse to a p-brane, 
arriving to the same equations of motion. 

Let S C X be the trajectory of a p-brane and let dS be the (p + l)-volume 
form of S. Rising the first index of dS and multiplying it by the tension t, we 
obtain a p-form with values on tangent vectors, that is called the impulse form 
of S. In other words, 

Definition 4.2. The impulse form of a p-brane S is the p-form on S with 
values on TS: 




U s : TSA p. ATS 



> TS C (TX)\ S 



defined by the following property: 



g(D , Us(D u D p )) = tdS(A), ...,D P ) 



for any Do, ... , D p tangent vectors to S. 
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If Do, . . . , Dp is an orthonormal frame of vector fields on S, where the matrix 
of g\s is diagonal (+1, -1, . . . , -1), then: 

p 

j=0 

where 5 = 1 and 5j = — 1 for j ^ 0. 

Example 4.3. In the case of a particle (p = 0), the trajectory S is a curve 
and the impulse form lis is a vector-valued 0-form; that is, it is simply a tangent 
vector lis = rnll, where m is the mass of the particle and U is the unitary tangent 
vector to the curve. 

Example 4.4. If (X = M} +n ,g = dt 2 — dxf) is the Minkowski spacetime, 
then the impulse form of a p-brane S can be written as: 

lis = ® dt + wi <g> d Xl + . . . + u n ® ^ n 

for some ordinary differential p-forms Ui on 5*. 

If S'tg := S H {t = t } is the particle at the instant t , then the vector: 

/ lis := ( / wo J 5 t + ...+ ( / u n )d Xn 
Js t0 \Js t0 J \Js t0 J 

can be understood as the total energy- impulse vector of the p-brane at to. 

The differential p-form uq is called energy form of the brane respect to the 
chosen inertial frame (t, x\, . . . , x n ), and the integral: 

/ w 

Js t0 

is understood as the total energy of the brane at t . 

Moreover, if the p-brane is at apparent rest at to (that is, d t is tangent to S 
at the points t = to) then it is easy to check the total energy of the brane at to is 
equal to: 

/ u — i - (Volume of S to ) ■ 
Js t0 

Definition 4.5. Let S G X be the trajectory of a p-brane, and let us write V for 
the Levi-Civita connection of (X, g). For any pair of tangent vector fields D,D' 
on S, the covariant derivative VdD' decomposes as a tangent vector to S plus a 
vector orthogonal to S: 

V D D' = tang(V DJ D') + nor(V D D') . 
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The first addend Vd-D' := tang(VD-D') is precisely the covariant derivative 
with respect to the Levi-Civita connection V of the submanifold (S,g\s)- 

The second addend is, by definition, the second fundamental form of S, 
which is a symmetric tensor with values on the normal bundle of S: 

$ s :TSxTS — > (TS)- 1 , $ S (D,D') := nor(V D D') . 

Therefore, the trace of the second fundamental form, tr $5, is a field of normal 
vectors to S. 



Proposition 4.6. The impulse form IT5 of a p-brane S satisfies: 

d v n 5 = dS <g> t-tr($ 5 ) . 

Proof. Let (D , . . . , D n ) be an orthonormal basis of vector fields on X, such that 
(D , . . . , D p ) is an orthonormal basis of vector fields on S. 
Let us write 



i=0 

where the Uij are the connection 1-forms. 



Consequently, 

n p 

V Dj Dj = ^ujijtDj) Di , and V Dj Dj = '^uj ij (D j ) D t , j < p 

i=0 i=0 

/ P \ P 

divs_Dj = contr (d^Dj) = contr J ^Ujj (g) Di J = u)ij(Di) , 



. i=0 / «=o 

where contr denotes the contraction of the contravariant and covariant indexes. 
Using these formulae, and taking t = 1: 

p p 

d v n 5 = d (^ d5 ) ® <W + (-!) P E(^d^) A 5,-dv^ 

i=o j=o 

p p n 

= ^(divsDj) dS ® ^-Dj + J2( iD i dS "> A ® Di 

j=0 j=0 i=0 

p p n 

i=o j=o i=o 
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(applying the formulae for div gDj and V DjDj) 

pp p 
= dS <g> Yl ,S j x 'i {I) ' ) ' I) j + Yl dS ® 5 ^ D i D i 

3=0 i=0 j=0 

P V P 

= -dS ®YY1 ^ji( D i) D j + dS®Y Sj^DjDj 

3=0 i=0 j=0 
P P 

= -dS g> ^ 8^7 Di Di + dS®Y SjVd^Dj 
i=o j=o 
p p 

= dS <g> ^ ^ ( V ^j " ^ A) = d5 ® 5^ fij$s(Dj, Dj) = dS ® tr $ 5 

5=0 j=0 



□ 



Example 4.7. In case p = 0, let 5 1 be the trajectory of a particle with impulse 
lis = mt/, where m is the mass of the particle and U is the future-pointing 
unitary tangent vector of the curve S. 
Since: 

$s(U, U) = nox{VuU) = VuU = V 9t U 
we observe tr $5 = Va r £7 is the acceleration vector of the particle. 

Definition 4.8. By analogy with the particle case just explained, if 5* is the 
trajectory of a p-brane, then the normal vector tr $5 is interpreted as the accel- 
eration of the brane. 

If there are no external forces, the trajectory S of a brane should have null 
acceleration. For a particle, this amounts to saying that it is a geodesic: Va T £7 = 
0. For a p-brane, this amounts to the equation: 



Inertial Motion: 



tr $c = 



By 14.61 this equation is equivalent to dyll^ = 0, which is an infinitesimal 
conservation law for the impulse. 

Remark 4.9. The equation tr$s = is precisely the Euler-Lagrange equation for 
the variational problem defined by the Nambu-Goto action. 
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4.2 Electromagnetic field 

Definition 4.10. An electromagnetic field over the spacetime X is a skew- 
symmetric tensor: 

F : TXA . P±\ . ATX — ► TX 
satisfying the following property: 

F(D ,...,D p ) e <D ,...,D p > x 

for any collection Dq, . . . , D p of vector fields on X. 

The value F(D , . . . , D p ) x may be understood as the force at the point x that 
suffers a brane with (p + 1)- volume vector Dq A ... A D p and unitary charge (see 
the force law below). 

The definition of F amounts to saying that the tensor: 

F(D , D p+l ) := g(F(D , ...,D P ), D p+l ) 

is a (p + 2)-differential form on X, and we will say that F is the (p + 2)-form of 
the electromagnetic field. 

Force Law for a p-brane 

Let S be the trajectory of a p-brane with tension t and electric charge q. 

Definition 4.11. The charge-current vector of this brane is the only (p+ 1)- 
vector Js on S satisfying 

dS(Js) = q . 

If (Dq, . . . , D p ) is an oriented orthonormal frame of vector fields on S, then: 

J s = q D A . . . A D p . 

Let F be an electromagnetic force and assume that the p-brane S does not 
substantially modify the electromagnetic field. Nevertheless, the p-brane S does 
suffer an acceleration due to the electromagnetic force F, that we postulate to be 
governed by the following equation: 



Lorentz Force Law: 



d v n 5 = ds ® f( j s 



Using Proposition I4.6[ this equation is equivalent to t ■ tr $ = F(J S ), which, 
substituting the value of Js, is in turn equivalent to: 

t-tr$ = q- F(D ,...,D P ) . 



20 



Observe the typical form of this equation: mass x acceleration = force; 
the definition 14.101 of F has been dictated by the need of giving sense to this 
expression. 

Example 4.12. In the case p — 0, the charge-current vector of a particle is 
simply a vector J$ = q U, where U is the future-pointing unitary tangent vector 
of the trajectory S of the particle. 

Since the impulse of the particle is mil, the force law reads: 

d v (mU) = dr®qF(U) 

where r stands for the proper time of the curve. 

As d^U = dr ® Vg r U, this force law is equivalent to the equation: 

m-V dr U = q-F(U) 

which is precisely the classical Lorentz Force Law for a particle of mass m and 
charge q. 

Remark 4.13. Let (X = ~R n+1 ,g = dt 2 — Y^i^x 2 ) be the Minkowski space-time. 
The trajectory of a p-brane S can be written as Xi = fi(t, u\, . . . , u p ), where 
(t, ui, . . . , Up) are local coordinates on S. 

On these coordinates (t,Ui, . . . ,u p ), the force law t • tr$ = F(Js) produces 
a system of second order partial differential equations, that is quasi-linear and 
hyperbolic. 

For these kind of systems, the Cauchy problem has a unique local solution 
([IS]. Proposition 3.2), so the force law uniquely determines the trajectory of the 
p-brane, for adequate initial conditions. 



4.3 Maxwell equations 

Definition 4.14. A distribution of charged p-branes on the spacetime X will 
be represented by means of a differential g-form C, called the charge-density 
form. 

The physical meaning of this g-form is the following (recall q = n — p): Given 
q linearly independent vectors Di,...,D q £ T X X, that we understand as an 
oriented infinitesimal parallelepiped at the point x, we have: 

C(D D ) — < ^ um ' ar I ec t e d with a sign, of the charges of the p-branes 
' ' ' ' ' q > y transversally crossing the parallelepiped 

We say that a p-brane with a trajectory S transversally crosses the paral- 
lelepiped Di, . . . , D q whenever T X X = T X S © (Di, . . . , D q ). If the orientation 
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of T X X coincides with the product of the orientations on T X S and (D\, . . . , D q ), 
then the charge of the p-brane counts with positive sign; otherwise, we affect the 
charge with a negative sign. 

Definition 4.15. The charge-current (p+l)-vector of a distribution of charged 
p-branes is the only (p + l)-vector J satisfying: 

ijdX = C . 

Equivalently, if J* is the (p + l)-form metrically equivalent to J and * stands 
for the Hodge operator, 

J* = (-l) pn * C . 

Example 4.16. When p — 0, the charge-density form C is a differential n-form, 
and the charge-current vector J is simply a vector field on X. 

In this case, the electromagnetic field F is a 2-form, related to the distribution 
of charges by the Maxwell equations: 

dF = , OF = J* . 



Let us consider a distribution of charged p-branes, represented by a charge- 
density g-form C or, equivalently, by a charge-current (p + l)-vector J. Such 
a distribution of charges "produces" an electromagnetic field, represented by a 
(p + 2)-form F. By analogy with the particle case, we postulate that both fields 
are related by the following: 



Maxwell Equations: 



or equivalently 



dF = 



dF = 



dF = J* 



d(*F) 



-lYC 



Remark 4.17. The second Maxwell equation implies an infinitesimal charge con- 
servation law dC = (equivalently, d J* = 0). 

Remark 4.18. The operators d and d are, essentially, the only first-order natural 
linear differential operators between differential forms. Therefore, in a certain 
sense, Maxwell equations are the only possible first-order equations that may 
arise. 

Variational principles. In a similar vein to what is done for charged particles 
(p = 0), the Lorentz force law and the Maxwell equations may be derived from 
variational principles, as follows. 
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Let us write F = dA where A is a (p + l)-form on spacetime, called the 
electromagnetic potential. For each trajectory S of a p-brane with tension t 
and electric charge q, consider the action: 

A(S) := -t [ dS + (-l) p q / A . 
Js Js 

Extremals of this action are precisely the trajectories that satisfy the Lorentz 
force law. 

On the other hand, let us fix a closed g-form C on the spacetime X. For any 
{jp + l)-form A, consider the action: 

A(A) := [ \-F A*F — [ A AC , 
Jx 2 Jx 

where F := dA. The Euler-Lagrange equations for this action amount to the 
Maxwell equation dF = J*, where ijdX = C . 

5 Fluid of charged £>-branes 

Now we shall extend in a natural way the notion of impulse of a p-brane and the 
force law to the case of a fluid of charged p-branes. Via ll.4[ we shall define the 
electromagnetic energy tensor associated to the electromagnetic field strength F, 
which is necessary to formulate the Einstein equation. 

5.1 Impulse form and force law for a fluid 

The following lemma is easy to check (v.gr. [11], Lemma 2.4): 

Lemma 5.1. The following linear map is an isomorphism: 

TX®TX — > K n X ®TX , T 2 H- C{ (dX <g> T 2 ) 

where C\ denotes the contraction between the first covariant and first contravari- 
ant indices. 

Moreover, ifT 2 is a 2-contravariant tensor on X and U n := C\{dX <g) T 2 ) is 
the corresponding vector-valued n-form, then: 

d v IT n = dX (8) divT 2 . 

If T 2 is a 2-contravariant tensor, we write T 2 for the 2-covariant tensor met- 
rically equivalent to it. 
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Definition 5.2. The mass-energy-momentum distribution of a fluid of charged 
p-branes is represented by a differential n-form IT n with values on TX, that we 
call impulse form of the fluid. 

The 2-covariant tensor T 2 corresponding to U n via the isomorphism of Lemma 
15.11 is called the stress-energy tensor of the fluid. 

The interpretation of the impulse n-form IT n is the following: assume the 
ambient manifold X is the Minkowski spacetime and let H be an oriented hyper- 
surface. 

If S is the trajectory of a p-brane transversally crossing the hypersurface 
H, and Us is the vector-valued impulse p-form of the p-brane, then the vector 
J SnH lis is said to be the total impulse of the particle in the hypersurface H (see 
Example EOJ • 

Now, the vector 



J H 

is understood as the sum of total impulses of all the charged p-branes transversally 
crossing the hypersurface H . 

Let us consider a fluid of charged p-branes, with impulse form Il n and charge- 
current (p + l)-vector J. 

In absence of external forces, the variation of the fluid impulse should be null: 
d v II n = 0. If an electromagnetic field is present we postulate, by analogy with 
the case of a single p-brane, that the movement of the fluid satisfies the Force 
Law: 



In virtue of Lemma [5.11 this equation is equivalent to divT 2 = F(J), or to 
divT2 = ijF . Combining it with the Maxwell equation dF = J*, we obtain 
another equivalent formulation: 



5.2 Example: Dust of charged j9-branes 

Let us consider a fluid of charged p-branes without pressure and where all the 
p-branes in the fluid have the same tension t and the same electric charge q. 
The general idea is that each p-brane has approximately the same velocity as the 
surrounding ones; hence, we give the following definition: 

A dust of p-branes is described by an integrable distribution on X of rank 
p+ 1, for which each integrable submanifold represents the mean trajectory of an 
infinitesimal portion of p-branes. 




d v II n = 



dX ® F(J) . 



divT 2 



i dF F 
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Let (Dq, . . . , Dp) be an orthonormal basis (+, —,...,—) of the distribution. 
Such a basis defines an orientation on each integral submanifold. By analogy 
with the case of a single p-brane, the charge- current (p + l)-vector of the dust is 
defined as: 

J := p e D A ... A D p , 

for some charge density function p e . 

The contravariant stress-energy tensor of the dust is defined by the formula: 

v 

T 2 := Pm J2 5 i D i® D i 

3=0 

where Sq = 1, 5j-^o = —1, and the function p m := (t/q)p e is called the mass 
density (that is, on each trajectory S of the dust we consider the dual metric 
(g\s)* multiplied by the function p m ). 

According 15.11 the corresponding impulse form is: 

v 

n n = p m ^2(i Dj dX) ® SjDj . 

3=0 



Proposition 5.3. If the charge conservation law dC = holds, then the impulse 
form of a charged dust satisfies: 

d v II n = p m dX®tr$ 

or, equivalently , 

divT 2 = p m tr$ , 
where $ is the second fundamental form of the trajectories of the dust. 

Proof. Let us complete the orthonormal basis (Do, . . . , D p ) of the distribution up 
to an oriented orthonormal basis (Do, ■ ■ ■ , D n ) of tangent fields on X, and let 
(9 , . . . ,0 n ) be the corresponding dual basis of 1-forms. 
The charge-density g-form is: 

C = ijdX = ij(9 A • • • A 9 n ) = Pe 6 p+l A • • ■ A 9 n . 

We have: 

v v 
n n = p m ^(i D AX) <g> SjDj = p m Yl ( e o A • ■ • A 9 P ) A (9 p+1 A • • • A 6 n ) <8> SjDj 

3=0 3=0 

= lD > (^o A • • • A ® 5jD 3 J A (p m 9 p+1 A • ■ ■ A 9 n ) . 
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and therefore, 



d v n n = |^d v ^ i Dj (0 O A ■ ■ ■ A 9 p ) <g> SjDjj A (p m 9 P+1 A ■ • • A 0„) 

+ (~ 1 ) P A ' ' ' A 9 p) ® S i D i) A d (P™ d P+i A • • • A 9 n ) . 



The second addend is null because d(p m 6 p+ \ A • ■ ■ A 9 n ) = dQC) = 0. With 
respect to the first one, the term which is differentiated has the same expression 
than the impulse form of each integral submanifold (considered as the trajectory 
S of a p-brane of tension 1). Applying Proposition 14.61 we obtain: 

d v IT n = ^d v { Oi (0o A ■ • • A B p ) ® 5jD^j A {p m 9 p+1 A • • • A 6 n ) 

= (9 A • • • A 9 p (8) tr $) A (p m 6 P+1 A • ■ • A 9 n ) = p m dX ® tr $ . 

□ 

As a consequence, if the electromagnetic field F satisfies the Maxwell equa- 
tions (so, in particular, the charge conservation law holds), then, for a dust, the 
force law divT 2 = F(J) is equivalent to: 

p m tr$ = F(J) . 
5.3 Electromagnetic energy tensor 

Definition 5.4. Let F be an electromagnetic field. Its electromagnetic energy 
tensor is the 2-covariant tensor T e \ m associated to the (p + 2)-differential form F 
according to Definition 11.41 

Let F be the electromagnetic field produced by a fluid of charged p-branes with 
stress-energy tensor T m . The Lorentz force law divT m = iqfF and Proposition 
11.71 produce an infinitesimal conservation law: 

div (T m + T clm ) = i dF F + (-i 9F F) = . 

Indeed, this property is the main motivation for the definition of the electro- 
magnetic energy tensor (see the Introduction). 

Finally, as in the particle case, we postulate that the electromagnetic energy 
has a gravitational effect through the Einstein equation. 
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To sum up, a fluid of charged p-branes is described by four tensor fields 
on spacetime: A stress-energy tensor T m and a charge-current (p + 1)- vector 
J representing the distributions of mass and charge, and a differential (p + 2)- 
form F and its energy tensor T e i m , representing the electromagnetic field and its 
electromagnetic energy. 

They are related by the following equations: 

Maxwell equations: dF = , dF = J* 

r(g) 

Einstein equation: Ricci(g) — g = T m + T c \ m 

Acknowledgements 

The authors acknowledge the referee for pointing out reference [5]. 



References 

[1] Anderson, I.M.: On the characterization of energy-momentum tensors, 
Gen. Relativity Gravitation, 10, 461-466, (1979). 

[2] Castrillon, ML, Munoz, J.: Gauge-invariant characterization of Yang- 
Mills-Higgs lagrangians, Ann. Henri Poincare, 8, 203-217 (2007). 

[3] Epstein,, D.B.A.: Natural tensors on riemannian manifolds, J. Differential 
Geom. 10, 631-645 (1975). 

[4] Gilkey, P. B.: Curvature and the eigenvalues of the Laplacian for elliptic 
complexes, Adv. Math. 10, 344-382 (1973). 

[5] Henneaux, M., Teitelboim, C.: p-Form Electrodynamics, Found. Phys., 
16, 7, 593-617, (1986). 

[6] Kerrighan, D.B.: Arbitrary Tensor Concomitants of a Bivector and a 
Metric in a Space-Time Manifold, Gen. Relativity Gravitation 13, 19-27, 
(1981). 

[7] KERRIGHAN, D.B.: On the uniqueness of the energy-momentum tensor for 
electromagnetism, J. Math. Phys. 23 (10), 1979-1980 (1982). 

[8] Kolar, I., Michor, P.W., Slovak, J.: Natural Operations in Differen- 
tial Geometry, Springer- Verlag, Berlin, (1993). 



27 



Lovelock, D.: The uniqueness of the Einstein-Maxwell field equations, 
Gen. Relativity Gravitation, 5, 399-408, (1974). 

Lovelock, D.: Bivector field theories, divergence-free vectors and the 
Einstein-Maxwell filed equations, J. Math. Phys., 18, 1491-1498, (1977). 

Navarro, A., Navarro, J.: Lovelock's theorem revisited, J. Geom. Phys., 
61, 1950-1956 (2011). 

NAVARRO, J., SANCHO, J. B.: On the naturalness of Einstein equation, J. 
Geom. Phys., 57, 1007-1014 (2008). 

Senovilla, J.M.M.: Superenergy tensors, Classical Quantum Gravity, 17, 
2799-2841, (2000). 

Slovak, J.: Peetre theorem for nonlinear operators, Ann. Global Anal. 
Geom., 6, 273-283, (1988). 

SLOVAK, J.: On invariant operations on a manifold with connection or 
metric, J. Differential Geom., 36, 633-650, (1992). 

SLOVAK, J.: On invariant operations on pseudo-riemannian manifolds, 
Comment. Math. Univ. Carolm., 33, 2, 269-276, (1992). 

Stredder, P.: Natural differential operators on riemannian manifolds and 
representations of the orthogonal and special orthogonal groups, J. Differ- 
ential Geom., 10, 647-660, (1975). 

Taylor, M.E.: Partial Differential Equations III, Nonlinear equations, 
Springer (1997). 



28 



